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TECHNICAL MEMORANDUM 1308 

ON MOTION OF FLUID IN BOUNDARY LAYER NEAR LINE 
OF INTERSECTION OF TWO PLANES* 

By L. G-. Loitsianskii and V. P. Bolshakov 


SUMMARY 

In the paper "The Mutual Interference of Boundary Layers , " the 
authors investigated the problem of the interference of two planes 
intersecting at- right angles on the boundary layers formed by the motion 
of fluid along the line of intersection of these planes. 

In the present paper, the results of the preceding one are general- 
ized to the case of planes intersecting at any angle. The motion of a 
fluid in an angle less than 180° is discussed and the enlargement of the 
boundary layers near the line of intersection of the planes, the limits 
of the interference effects of the boundary layers, and the corrections 
on the drag are determined. All computations are conducted by the 
Karman-Pohlhausen method for laminar and turbulent boundary layers. The 
results are reduced to tabulated form. 


INTRODUCTION 

The problem of the Interaction of the boundary layers formed in the 
dihedral angle between two thin plates parallel to the intersection of 
the plates was apparently first proposed in reference 1 (for the case 
of a right angle). In the discussion of reference 1, the need for a 
generalization of results obtained fpr the case of the right angle to 
the case of any angle less than two right angles (180°) was pointed out. 
The present paper is concerned with this problem and its solution. 
Although the limits of application for. the approximate method of the 
finite -thickness layer previously used are retained, the problem of the 
interaction of the boundary layers near the intersection of a dihedral 
angle of any magnitude from 0° to 180° is solved herein. For a laminar 
layer, a first, and a second approximation are given, and also, for a 
check, a sixth approximation ( in the terminology of Pohlhausen). It is 

/ 

*"0 Dvizhenii Zhidkosti v Pogranichnom Sl6e Vblizi Linii Peresechenia 
Dvukh Ploskostei." Rep. No. .279, CAHI, 1936, pp. 3-18. 
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shown that the sixth approximation differs comparatively little from the 
second. In conclusion, the case of a turbulent boundary layer is considered 
with the assumption of the validity of the ' 1 / 1 ' power law for the velocity 
profiles . As in reference 1, all computations can be carried through to 
the end, although the procedure is somewhat more cumbersome. The limits 
of interference of the layers and the correction on the drag due to the 
interference effect are determined. 


1. DERIVATION OF FUNDAMENTAL INTEGRAL CONDITION 

Consider the flow of a fluid, approaching from infinity with 
velocity V, in the dihedral oblique angle 6 between two plates of the 
same finite length x along the flow and infinite in the transverse 
direction (fig. l) . The Y- and Z-axes are taken along the leading edges 
of the plates and the X-axis along the flow parallel to the line of 
intersection of the plates . In the oblique system of coordinates thus 
obtained, the distribution of the velocities in the boundary layer may 
be given in the same manner as for the case of the rectangular system 
of coordinates for the flow in a right angle. 

If each plate worked independently of the other, there would be 
formed on it a layer of thickness 8q, which is a function only of x 
and the profile of the longitudinal velocities 

u 0 = z > 3in 


or 


u 0 = u o( x > ain 0 ) 

depending on whether the boundary layer is considered to lie in the 
plane XOY or XOZ. 

Because of the retarding effect of one of the plates on the other 
near the line of intersection, the layers must not be considered as in 
the plane problem. The layers become three dimensional and the thick- 
ness will now be a function of two variables : 

6 »= 5(x,y) on the plate containing the Y-axis 

5 = 8(x,z) on the plate containing the Z-axis 

The component of velocity parallel to the X-axis, which will be denoted 
by u without the subscript 0, will be a function of the three 
variables ; that is , 


u = u(x,y,z) 
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By the conditions of the problem, under the basic assumption of the 
concept of a finite region of influence of the viscosity, this function 
must become zero on the surface of the plates, a constant at the outer 
limit of the boundary layer, and, particularly, at a finite distance from 
the line of intersection of the plates, must become the velocity distri- 
bution uo(x,z, sin 0) or uo(x,y, sin 0), which corresponds to the 
isolated plate with boundary layer undisturbed by the adjacent plate. 

The boundary separating the region disturbed by the adjacent plate from 
the undisturbed region, which corresponds to the isolated plate, shall, 
for briefness, be denoted as the "interference boundary" of the layers. 
The equation of this boundary in the planes XOY and XOZ will be 

h 0 = h 0 (x) 


A section of the boundary layer cut by the plane x = £ is shown 
in figure 2 . Inasmuch as the coordinates of the system YOZ are 
oblique, the equation of the boundary layer in this section for the 
undisturbed region for y>y 2 (5)> where 

y 2 (€) = 003 9 = V e) ” 5 0 (€) Got 9 ’ wil1 have thB form 


z „= z iU) = 


& 0 (O 

sin 0 

\ 


Equation (l) holds true both for 0<g- anh for 



( 1 ) 


The non- 


coincidence for 0 / of the interference boundary h Q (£) an ^ 

coordinate 7q(%) should be noted, from which starting point the 
ordinate z of the outer limit of the boundary layer becomes and remains 
constant, independent of y. All that has been stated about the YOX 
planes also remains true, of course, for the ZOX plane (because of 
the symmetry y^(g) = and y 2 (g) = z %(K) > which should be remembered 

in the following development) . 


For the present, the question as to the equation of the boundary 
layer in the section x = will be ignored. The fundamental integral - 
condition of the problem will now be set up. For this purpose, as in 
the work previously cited, the momentum theorem is applied for a tube 
of flow formed by the coordinate planes, by the surface of the lines of 
flow passing through the edge of the boundary layer at section x *= x, 
and, finally, by the surfaces of the lines of flow passing through the 
perpendiculars to the plates located at the points y = hg(x) and 

z = hg(x), where the part of the flow tube considered is between the 

sections x = x and the section located upstream of the flow at a 
sufficient distance from the section x = 0, that is, from the entry of 
the fluid on the plates . Then, as is known. 
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p(V-u) u da = W 


( 2 ) 


where W is the drag of the plates applied to the segment of the flow 
tube considered and a is the section of the boundary layer cut by a 
plane perpendicular to the X-axis at distance x from the origin 0. 


The double integral on the left may be expressed in the following 
manner (fig. 2): 




■where Z, the ordinate of the outer edge of the boundary layer, is, as 
yet, an unknown function of y and x if the angle 6 is considered 
as a par ame ter maintaining a constant value for the given problem. The 
first two integrals have a very obvious originj the presence of the last 
integral is due to the obliqueness of the coordinates and the differences 
in the direction between the ordinate zp(x) and the thickness of the 
layer Sq(x) , which makes it necessary to take the integration over the 
two areas of the triangles shown hatched in figure 2. 

The drag W, as is easily seen from figure 3,- is determined by the 
integral 
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Lvio 


•I 


h(x) 


T dy + \ T„ dy 

V« 


n x 


d£ * 2 


P72(5)+Zi(€) coa 0 




^y 2 (x)+z 1 (x) cos 0 


T n 


Jy2(^) +z i(C) coa 0 


d£ 


'4) 


where T denotes the friction stress in the region S disturbed by the 
adjacent plate and denotes the corresponding stress in the undis- 

turbed region Sq . The boundary between the regions S and Sq is , 
of course, the boundary of interference of the boundary layers . 

Combining equations (3) and (4) yields, in general form, the 
integral condition of the problem: 


sin 0 


= 2 


x 


vJo 


r 

•i 

i 


v 2 


r z 


dy 


p( V-u) u dz + 


dy 


p(Y-u) u dz + 


y 2 +z^ cos 0 


y 2 +z-|_ cos 0-y 


cos 0 


dy 


p(V-u) u dz 


y 2 (C)+z 1 (?) cos 0 


I 


y 2 (x)+z 1 (x) COS 0 


t dy + I dy 

)y 2 (C)+ z i(?) cos 0 


d£ 


(5) 


Equation (5) is a generalization of the integral condition derived in 
the previously cited paper for the case of a right-angled dihedral 

(0 = —)' '. The limits of integration are as yet unknown functions of the 

arguments. In order to determine the form of these functions, it is 
first necessary that the shape of the boundary layer be given. 


T dy + 
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2. SHAPE OF CROSS SECTION OF BOUNDARY LAYER; 'VELOCITY DIAGRAMS; 

SOLUTION OF PROBLEM BY FIRST APPROXIMATION 

The section of the boundary layer on the plane x = £ (fig* 2) is 
considered and it is assumed that the equation of the curve defining the 
edge of the "boundary layer is 


f(y,z) = a(£,0) 

y 1 (4) - i - y 2 (Z) 

Z 1 (€) - z - z 2 (C) 

The form of the function f(y,z) cannot be determined unless the 
additional assumption is made as to the similarity of the approximate 
velocity diagrams in the different planes x = The curve is sought 
in the form 



where the function <p(t) is subjected to the conditions 




and the function f(y,z) is subjected to the conditions 

f(o,z) = f(y,o) = 0 (9) 

Then the diagram of velocities of equation (7) will evidently satisfy 
the end conditions of the problem at the walls and at the edge of the 
boundary layer of equation (6). Setting 

y = iU,e) y' 
z = i(z,e) z' 

where 2(£,©) is a certain length characteristic for the section x = £ 
and y' and z' are nondimens ional magnitudes, yields equation (7) in 
the form 

u 

y = <P 
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From the condition of similarity of the velocity diagrams, the right side 
of this equation must not contain £, which can he the case if f(y,z) 
is -a homogeneous function of the n^h degree, so that . 

where 


1 — = constant 

a(5;0) 


( 10 ) 


the constant depending, of course, on 0. It is easy to see that the 
degree of homogeneity should he equal to 2 (n = 2) and the function 
f(y,z) must simply he equal to the product of the variables yz because 
otherwise, with conditions (9) satisfied, the derivatives of the velocity 
with respect to the coordinates would become zero at the walls and would 
not give any friction. From equation (6), the following equation of the 
edge of the boundary layer at the section x = £ is obtained: 

yz = a(5,0) (11) 

where, from equation (10), for example, choosing the coordinate 
y 1 (g) = zq(£) of the undisturbed layer for the characteristic length 

2(£,0) yields 


a(e,e) = k(e) z 2 (5) = -M|L 8 2 (S) (12) 

x sin^ 0 

where k(0) is a certain nondimens ional constant parametrically 
dependent on the angle 0. 

It is now easy to obtain the magnitudes y 2 (£), 1 iq(£), and the 
function Z (y,?,0) for given angle 0. From equations (ll) and (12), 

■*!<«> - t^ftf ‘ k(&> ‘ 1<V * Ml 2 ? 8 0(O (13) 

where 5 q(£) is a known function of £ independent of the angle 0 
and determined by solving the problem of the boundary layer of the 
isolated plate. 

From the definition of the Interference boundary ho(g) given in 
equation (l), 
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h 0 (Z,e) = y 2 («) + z x (€) cos & = [k(S) + cos ejzjte) = - (0 3 ^ - g Oa 9 5 0 (g) 

(14) 

Finally, from equation (ll), there directly follows that 


z(y,g,e) 


a(g,Q) _ k(0) _ Sp 2 U) 


y 3 in^ 0 


(15) 


It is easily verified that, for y = y 2 (£), Z "becomes z 1 (g). The 
velocity diagram in the disturbed region of the "boundary layer will 
therefore "be smoothly Joined with the velocity diagram in the undisturbed 
region, that is, with the diagram of the isolated plate. 


When the equation of the "boundary layer ia found, the required 
velocity profile is obtained from equation ( 7 ) ; 



or from equation (15) : 



If Y denotes the variable ordinate of the edge of the boundary layer, 
that is, the magnitude that, from equation (ll), satisfies the equation 


Y-Z = a(£,e) 


then 



(17') 


It is important to note that equations (16), (17), and (17') are true 
not only for those values of y and z that satisfy inequalities (6): 

y 1 (5,0) - y- y 2 (5,e) 


z 1 (?,0)- z »z g (C,0) 


but in the entire range of interest; 
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oSy«y 2 (£,0) 

os Z Sz 2 (£,e) 

When the velocities of the points located in the rhomb 

0<z< Zi(5,0) . 

are considered, the boundary layer for these points is as though infinite, 
but the velocities are determined from computation on the hyperbolic edge 
of the boundary layer. 

The velocity profile has thus been determined and the edge of the 
boundary layer is known. Substituting the values of u and Z in the 
integral condition (5) yields an ordinary differential equation with one 
unknown a(£,0) inasmuch as all the magnitudes involved, including the 
friction at the wall 





are expressed in terms of this function. 

As is easily seen, however, the differential equation reduces to 
a simple equation in finite form for determining the coefficient k(0) 
appearing in equation (12). In order to obtain this single unknown 
coefficient, certain boundary conditions are assigned for the function 
<P(t) and Its derivatives, as in the classical Karman-Pohlhausen method. 

A consideration of the first approximation is the first step. The 
function <p(t) is subjected only to the conditions 

q»(6) = o 

<p(i) = i 

that is, the velocity u(x,y,z) becomes zero at the walls and 1 at the 
outer limit of the boundary layer, which leads to the profile 
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u _ yz 
V a(x,0) 


(19) 


For substitution in the integral condition, the integrals that enter 
this condition are first computed: 



z dz = 


V 

a(x,0) 


y 2«j 2 

- T a(x,0) Y 


r. 

Jo 


py2 /-> z 

dy | u dz s= 

Jyp 


0^2 

dy 

oy x vJ 


T r 


r 

0 

^2 


u dz 


a 2 (x,0) 


a (x,0) .| yi 2y 


dy « |-Va(x,e) logg 

J 1 


r\y 2 + z i 003 9 


y 2 + z i cos ©-y 


cos © 




u dz 


O h n(x) P 


= 2 


V 


i(x,0) 


"0 

y dy 

y2 vjo 


h 0 (x)-y 

cos 0 


z dz 


r\ h 0 ( x ) 


aCx-^©) cos 2 ©J 
a(x,0) jj a ^ X > 


^2 


[* 


h 0 2 (x)-2h 0 (x) y + y z y dy 


0) cos 0 + Yg y 


cos 2 ©J 
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In a similar manner, the integrals of the square of the velocity are 
computed. The development will he limited to computing the integrals 
just given and the remaining ones are written out in ready form: ■ 


‘71 


t 


dy 1 u 2 dz = ~ V 2 a(x,0) 

vJo Oo y ' 


I * y 2 

dy 

yp . lo 


a (xp) 0 


u 2 dz = ~ V 2 a(x,0) log 0 


,h n (x) 


, h o( x )-y 

cos 0 


dy 


* 


.2 1 2 


i yp 4 cos 2 a 


u dz = Y [e y x cos e + 15 a(Xj0) + 
P y x 6 cos 3 0 "| 

a 2 (x,0) J 


Thus, the left side of the integral condition reduces to the form 




vi vi 


(cr) p ^ V_u ^ u da = pY 2 sin 0^ a(x,0) + i a(x,0) log e 


^ y^ 1 cos^ 0 ^ y^ b cos’ 5 0 

6 y l 2 003 0 + 60 a(x,0) 90 a 2^ 0 ) 


] 


The computation of the right side of integral condition (5) reduces to 
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n x 


OO 


] 

lOo 


h (g) ryh (x) 

0 V ' r 0 

T ay + l dy 

0h 0 (O 


o * 


dg = ~ 

3 a n 


sin 0 


U0 


I 


h 0 (S) 


[i V 




ciy 


V x) 


h>«> 


“ T 50?7 


^ 1 ^cos 2 e 

sin 0 1 


y x 2 (g) 

a(g,0) 


ag + 


^ + 003 8 1 5#-£ ^ 


tfy a5 


According to the first approximation for the Isolated plate (according 
to Pohlhausen ) , 


V«) 

»o«) 


*i<«> - n 


i-A/i 

In 0 V 


12 vg 


Ox 


vJ 


sin 0 sin 0 V T 

O x 


ag . / V 

jysT ” sln e V i2» , 

O’ 1 \J0 


ag 1 V ^ . . 

— = ? - S Q (x) sm 0 


so that 


W = 


_ J_X. 


sin 0 


r» x 


2 o 

cos 0 


\J0 


r> x 


yi (g) 

a(g,e) 


d g 


vJO 


a(g,e) 

y x 2 (g) 


dg + 


1 a(x ; 0 _ ) _ y 3ln 2 e + - S 0 2 (x) - cos 0 

O v O ^ V 


V 
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All these expressions are Immediately simplified as soon as equation (12) 
is considered; then 


a(5,e) = — 5 0 2 (5) 
sin^ 0 u 

After all the simplifications are made, the following equation for 
determining k( 0 ) is obtained: 


2 3 

i . /„\ 2 cos 0 2 cos 0 1 cos 0 

iog e k( 0 ) - 3 + + 5 + is 


k3(0) 


( 20 ) 


It is readily seen that, for 0 = £-, equation (20) becomes 


log k = £. 
e 3 


k = 1.95 


as given in the previously cited work on the interference of the boundary 
layers on mutually perpendicular plates . 

The explicit dependence of k on 0 is not given in equation (20). 
This relation may be obtained by. the following simple device: When 


cos 0 v 

k(0r ' * 


(21) 


equation ( 20 ) becomes 


cos 0 


£e 


— + y + — y 2 

3^5^ 


1 w3 
15 4 


( 22 ) 


When the values of £ are given in the interval where the absolute value 
of the right side of equation ( 22 ) will not exceed 1 , the corresponding 
value of 0 is obtained, Y, is determined, and then k( 0 ) is found 
from equation (21). It would also be possible, of course, to proceed in 
another manner: Equation (22) may be rewritten 


k(S) 


3 ^ 5 ^ 


+ -L w3 
15 


( 22 ') 


From a given value of Y, , k(£) is obtained and then cos 0, and so forth. 
The simplest method is to draw the function 
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2 + y , 2 2 _L 

3 ^ + 5 ^ + 15 ^ 

T) = *(£) « (23) 

and obtain its intersection with straight lines parallel to the £-axis 

T| - cos 0 


and then to obtain k from K- 

The values of 0°, k(0), and £(0) are tabulated in table 1. 

The dimensions of the region of interference of the layers are first 
determined. According to equation (14), 

h Q (x ? 0) = m(0) 8 q (x) (14') 

where 


m 


k + cos 0 1+K 

sin 0 = cot 0 


(24) 


The values of m for different 0 can be determined from table 1 
The value 0 = 180° is somewhat isolated; for this value of 0 the 
value of m becomes indeterminate. The value of m can easily be 
obtained, however, by the usual device of analysis : 



-fe 5 cos 2 

\tan 0/ 0=rt 

\d 8 


e) ..f 

' L? 


'C Y 2 (£) sin 0 
Y(fc) + £ ¥' (£) 


0=1 


0 


so that, for 0 = 180°, K = - 1 and m = 0. The values of m(©) are 
given in the last column of table 1. 


The correction in the drag of one side of the plate due to inter- 
ference of the layers may be computed by the equation 



r 1 

1 — 

w 

o 

A 

C__ 

AW = 


Jo (T °' T> dJ J 

V 

0 


(25) 


inasmuch as the difference in the drag may show up only in the region S 
(fig. 3). The quantity l is the length of the plate in the direction 
of the flow. 
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The integral equation (25) is readily computed and yields 

^(c) +z i (£) 008 6 


AW = 



(t q -t) dy 


cL£ 


Finally, 



AW = p(0) p Yl 


p(0) = Lil: COB 2 9 

2k(0) sin 0 J 


(26) 


If L denotes the width of the plate in the direction transverse 
to the flow and the effect of its free end is not considered, the 
relative correction due to the interference of the boundary layers may 
he computed : 


aw _ p(e) u yi p(0) z 1 nr’ 

¥ = 0.578 L °- 578 L V™ 

or, when the Eeynolds number of the plate, which is equal to Yl/v, 

is introduced, 


AW 

W 



where q(0) denotes the relation 


q(e) 


p(Q) 

0.578 


(27) 


Generally, the correction obtained is extremely Insignificant for plates 
that are long in the transverse direction. If, however, the transverse 
length L is comparable with the width of the region of the disturbed 
layer at the end of the plate x =2, the correction is not insignificant. 
Thus , if 
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L = n*liQ(Z) « n*m(0) 8 q(Z) = n*m(9) \J — A 

the relative correction will now he equal to 

AW _ q(0) 1 

¥ Vl2 m(9) n 


or 


where 


aw fs(e)~| P ercent 
W =Ln J 

„ 100 l(9) 

Vl2 m(0) 


(28) 

(29) 


All the magnitudes introduced in the preceding equations may be expressed 
in terms of the previously given parameter £ : 

v(K) = 

q(S) = 


s(S) = 

In table 1, it is possible to find the values of these magnitudes for 
different angles 0 or the corresponding values of the parameter <£. 

The magnitude s (0) increases with an increase of the angle 0. Some- 
what paradoxical is the value of 50 percent for s(0) at 0° = 180 
and m = 0. It is found that the region of interference is equal to 
zero, and that a relative effect occurs on the drag, which is due to the 
fact that the width of the region and the absolute correction on the 
drag simultdneously approach zero. On the other hand, when 0 decreases 
to 0 and m->®, s decreases to 16 percent. Both these cases are 
limiting cases. 


1-^ 

cot 0 

1—5^ 

1.156 K cot 9 ) 
1 -y 

■ 100 

1.156 Vl2 


(30) 


In figure 4, the curves of the relation between X> and 0 are 
drawn for the first and other approximations and also for the case of 
turbulent layers . 


■ ii i 
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3 . SECOND AND SIXTH APPROXIMATIONS 

The following approximations differ from the first only in the 
boundary conditions that are imposed on the function 9 (t). 

The second approximation corresponds to the assumptions 

9 ( 0 ) « 0 

9(1) = 1 

9 ' (l) - 0 


The velocity profile for these conditions will have the form. 


u = 



No new difficulties in principle, as compared with the first approximation, 
are obtained. Inasmuch as the basic computations have been explained in 
the first approximation, the results are presented with the same 
notations . 


The magnitudes 
The final form 

log e k(0) = iSj 


are all given In table 2. 

of the equation for the determination of k(©) Is 

3 cos 9 15 cos 2 6 2 cos 3 6 5 cos^ 9 

+ TW + 28 t 2 (e) + 2l‘ k 3 (e) - 168' t 4 (9) + 


1 cos 3 9 
42° 


(31) 


The solution Is effected by the same device as in the case of the first 
appr oximat ion. 


All computations were also carried out for the sixth approximation, 
by which Is meant the results of the assumption of the following 
boundary conditions: 


min inii ■ i m 
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9 ( 0 ) « 0 f(l') = 1 

*"(0) =0 <P'(1) m 0 

9 "'( 0 ) « 0 <p"(l) = 0 

<p'"(l) «= 0 

Tills approximation gives, for the undisturbed layer, the velocity 
profile 

‘o-’f fe) - 5 fej + - 2 ©1 

and the following values of the frictional force and layer thickness, 
which are extremely close to the accurate solution of Blasius: 

T 0 - 0 - 33 ^ 

s o * 

It is therefore reasonable to suppose that, for investigating the prob- 
lem of the interference of the layers, the sixth approximation is 
considerably more satisfactory than the second. 

The results are again collected in table 3. 

In figure 4, the curve of K against 9 is given for the second 
and sixth approximations . 

As before, the equation for the determination of k(0) is also 
given 




p pCI ^ Q 

0.0524 + 0.0053 

k 4 


cos^ 9 


+ 0.0025 


cos® 9 


( 32 ) 
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4 „ raTEEFERENCE OF TURBULENT BOUNDARY LAYERS 

For the solution of this concluding part of the problem, the well- 
known ’1 /l' power law is applied, which, when the boundary layer is 
not two dimensional, assumes the form 



As is known, for an infinite plate without interference of the 
layers , this law gives excellent agreement with experiment ; it may he 
expected that the application of the ' l/7' power law to this case will 
he justified hy experiment . 

The generalized integral condition (5) is substituted in the 
velocity-profile equation (33) . Unfortunately, the integrals will now 
not he so easily evaluated because the presence of fractional exponents, 
particularly in the integrals taken over the hatched areas of the 
triangles shown in figure 2, strongly complicates the computations and 
leads to the necessity of taking integrals of binomial differentials. 
All the computations can be made with a sufficient degree of accuracy, 
however, by using converging series and integrating them. The left 
side of the integral condition reduces to the following expression: 



In computing the right side of the integral condition, the 
Kara nan formula is used for the expression of the friction at the wall 
In the undisturbed region of the layer: 


T 0 = 0.0225 PT 2 ^ l/4 = 0.0225 pV 2 ^ “ ln ^ 


(35) 
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Tn the disturbed part of the layer, the analogous formula 
T. 0.0225 PV 2 ^ 4 - 0.0225 

■ °' 0225 cT3 (v a(x%! Bln e) lA < 36 > 


is assumed. These values are substituted, aa in the previous sections, 
in the formula for the resistance of the vails forming part of the 
boundary of a tube of flow. Then, 
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The Kerman formula was then used: 


® 0 M - 0 ' 370 fe) 1/5 


When It is known (see equation (12)) in advance that the solution of the 
problem will have the form 


2 ' 5 0 W 

a(x,0) - k(0) y (x) = k(0) 5 — 

sin 0 


(38) 


this value is substituted in equations (34) and (37); then, after some 
simplifications, the following solution is obtained; 


fl 


i 2 b 0 z (x)( 

p(Y-u) u do « p — — lo . 1607 k + 0.0972 k log Q k + 0.1360 cos 0 - 


sin 0 


0.0230 


W = p 


cos 2 0 


k 

^5 2 (x) 


+ 0.0023 


cos 3 0 


- 0.0004 


cos^ 0 


+ 0.0001 


cos 5 0 






'0 / non 2 0 

'0.1949 k + 0.1949 cos 0 - 0.0195 k + 0.0122 — — Z 

1- 


sin 0 


0.00305 


cos 3 0 


+ 0.00133 


cos 1 ^ 0 
k 3 


+ 0.00073 


cos 5 0 


+ • 


J 


Equating these two expressions yields the following equation for 
determining k(0) : 

log k = 0.1518 + 0.6056 C °y ° + 0.3632 Coa ~- 6 - 0.0556 C03 ^ — + 
e k yd ^3 


_ COS^ 0 _ orvoo COS 5 0 „ nnn COS 6 0 

0.0175 ; 0.0082 = — + 0.0051 = — 


+ • 


As in the preceding sections, the change of variables is made: 


cos 0 

kTeT 


= K(e) 


(39) 
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Tlie transcendental equation then becomes , 

cos 6 = X exp (0.1518 + 0.6056 X + 0.3632 X 2 - 0.0556 X 4 5 + 

0.0175 'd - 0.0082 X 5 + • • •) (40) 

This equation Is easily solved by the tabular method, where 

0°< 9< 180° 0.5540 >£>- 1.0000 

The corresponding values of Q, k(0), and £(0) are given in table 4. 
The further investigation in no way differs from that for laminar layers. 
The values of the coefficients characterizing the boundary of the region 
of interference and the corrections on the drag are given in table 4 
for various values of 0, with all the computations ommitted, in the 
notation previously used. 


In connection with the formulas of turbulent friction, the 
coefficients q and p are determined by the formulas: 


and 


AW *« p*p 



(41) 



(41*) 


The dependence of p on X is determined by the following series : 
p _ i C ot 0 (0.00133 - 0.00083 X 2 + 0.00021 X 5 - 0.00009 X 4 + 


0.00005 X^ - 0 . 00003 X ® + . . .) 


(42) 


4 ~ 0.036 


(43) 


The dependence of X on 9 is given for the case of turbulent layers 
in figure 4. The boundary of the region of interference (the 
coefficient m(0)) in the case of the turbulent layers differs little 
from the corresponding boundary of the la m inar layers, according to 
the sixth approximation; whereas the relative correction on the drag s 
in percent is several times less in the turbulent case. 
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CONCLUSION 

The results obtained have a readily understandable form.. In general, 
the effect of the interference of the layers on the drag of the plates 
is insignificant. The effect assumes an appreciable value only in the 
case where the plates in the dimensions transverse to the flow become 
comparable with the width of the region of disturbed boundary layer. 
Moreover, interference plays a large part in the motion of fluids through 
small dihedral angles. Thus, for example, in the motion near the inter- 
section of a dihedral angle of about 10°, the region of interference 
exceeds by 16 times the thickness of the layer at the given section. 

At smaller angles , the phenomenon is still more marked . 

All the conclusions of the present and preceding papers require 
experimental check. 

By agreement with the Central Aero -Bydrodynami cal Institute, the 
aerodynamic laboratory of the Leningrad Industrial Institute is under- 
taking an experimental investigation of the' phenomenon of the inter- 
ference of boundary layers . It is proposed, through use of bhe method 
of microtunnels, to observe directly the distortion in the velocity 
profiles, and so forth, of the phenomenon. 

The present work was carried out at the Aerodynamic Laboratory of 
the Leningrad Industrial Institute. 


Translated by S. Reiss, 
National Advisory Committee 
for Aeronautics 
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TABLE 1 


Q 

(deg) 

k(0) 

W) 

m(6>) 

P(e) 

1(0) 

s(S) 

(percent) 

0 

.2.894 

0.3455 

*» 

am 

m 

16.3 

10 

2.880 

.3419 

22.257 

7.324 

12.671 

16.4 

20 

2.840 

.3309 

11.315 

3.697 

6.396 

16.7 

30 

2.773 

.3123 

7.278 

2.503 

4.330 

17.2 

40 

2.682 

.2856 

5.364 

1.916 

3.315 

17.8 

50 

2.567 

.2504 

4.190 

1.570 

2.716 

18.7 

60 

2.434 

.2054 

3.388 

1.346 

2.329 

19.8 

70 

2.283 

.1498 

2.793 

1.187 

2.054 

21.2 

80 

2.120 

.0819 

2.329 

1.069 

1.849 

22.9 

90 

1.948 

0 

1.948 

.974 

1.685 

25.0 

100 

1.773 

-.0980 

1.624 

.892 

1.543 

27.4 

110 

1.601 

-.2136 

1.340 

.813 

1.407 

30.3 

120 

1.442 

-.3467 

1.088 

.732 

1.266 

33.6 

130 

1.300 

-.4945 

.858 

.641 

1.109 

37.3 

140 

1.185 

-.6464 

.6 2 

.536 

.927 

41.0 

150 

1.099 

-.7880 

.466 

.417 

.721 

44.7 

160 

1.042 

-.9018 

.301 

.285 

.493 

47.3 

170 

1.010 

-.9750 

.144 

.142 

.246 

49.3 

180 

1.000 

- 1.0000 

0 

0 

0 

50.0 


TABLE 2 


9 

(deg) 

k(e) 

K(e) 

m(0) 

P(0) 

l(0) 

s(e) 

(percent) 

0 

2.217 

0.4511 

m 

m 

m 

13.7 

10 

2.205 

.4466 

18.360 

10.162 

13.921 

13.8 

20 

2.168 

.4334 

9.090 

5.148 

7.052 

14.2 

30 

2.108 

.4108 

5.948 

3.584 

4.910 

14.7 

40 

2.026 

.3781 

4.343 

2.702 

3 .701 

15.6 

50 

1.925 

.3339 

3.352 

2.233 

3.059 

16.7 

60 

1.808 

.2765 

2.665 

1.928 

2.641 

18.1 

70 

1.678 

.2038 

2.150 

1.712 

2.345 

19.9 

80 

1.542 

.1126 

1.742 

1.546 

2.118 

22.2 

90 

1.407 

0 

1.407 

1.407 

1.927 

25.0 

100 

1.284 

-0.1353 

1.127 

1.279 

1.752 

28.4 

110 

1.182 

-0.2893 

.894 

1.153 

1.579 

32.2 

120 

1.108 

-0.4513 

.702 

1.018 

1.395 

36.3 

130 

1.060 

-0.6064 

.544 

.875 

1.199 

40.2 

140 

1.032 

-0.7422 

.414 

.721 

.988 

43.6 

150 

1.016 

-0.8524 

.300 

.556 

.762 

46.3 

160 

1.007 

-0.9332 

.196 

.379 

.519 

48.3 

170 

1.002 

-0.9832 

.097 

.192 

.263 

49.6 

180 

1.000 

- 1.0000 

0 

0 

0 

50.0 
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TABLE 3 


e 

(deg) 

k(e) 

Z(e) 

m(e) 

p(e) 

i(e) 

s(e) 

(percent) 

0 

1.969 

0.5079 

• - ■ 

m 

m 

12.3 

10 

1.957 

.5032 

16.940 

8.416 

12 . 713 

12.4 

20 

1.921 

.4892 

8.365 

4.272 

6.453 

12.8 

30 

1.863 

.4648 

5.458 

2.921 

4.412 

13.4 

. 40 

1.783 

.4296 

3.965 

2.262 

3.417 

14.3 

50 

1.685 

.3815 

3.039 

1.879 

2.838 

15.5 

60 

1.572 

.3181 

2.393 

1.631 

2.464 

17.0 

70 

1.450 

.2359 

1.907 

1.457 

2.201 

19.1 

80 

1.327 

.1309 

1.524 

1.324 

2.000 

21.7 

90 

1.204 

0 

1.204 

1.204 

1.819 

25.0 

100 

1.106 

-0.1571 

.946 

1.094 

1.653 

28.9 

110 

1.038 

-0.3295 

.741 

.985 

1.480 

33.2 

120 

1.004 

-0.4982 

.582 

.887 

1.340 

37.5 

130 

1.000 

-0.6428 

.466 

.766 

1.157 

41.1 

140 

1.000 

-0.7660 

.364 

.643 

.971 

44.2 

150 

1.000 

-0.8660 

.268 

.500 

.755 

46.7 

160 

1.000 

-0.9397 

.176 

.342 

.517 

48.5 

170 

1.000 

-0 . 9848 

.088 

.174 

.263 

49.6 

180 

1.000 

- 1.0000 

0 

0 

0 

50.0 


TABLE 4 


e 

(deg) 

k(9) 

Z(d) 

m(0) 

p(e) 

q(0) 

s(e) 

(percent) 

0 

1.805 

0.5540 

m 

«• 

• 

5.35 

10 

1.796 

.5483 

16.010 

0.0115 

0.3189 

5.38 

20 

1.767 

.5318 

7.915 

.0058 

.1614 

5.51 

30 

1.720 

.5035 

5.172 

.0039 

.1094 

5.73 

40 

1.654 

.4631 

3.765 

.0030 

.0839 

6.03 

50 

1.577 

.4076 

2.898 

.0025 

.0689 

6.43 

60 

1.485 

.3367 

2.292 

.0021 

.0594 

7.00 

70 

1.382 

.2475 

1.835 

.0019 

.0525 

7.73 

80 

1.273 

.1364 

1.469 

.0017 

.0472 

8.68 

90 

1.164 

0 

1.164 

.0015 

.0431 

10.00 

100 

1.065 

-0.1631 

.905 

. 0014 

.0394 

11.76 

110 

1.000 

-0.3420 

.700 

. 0013 

.0364 

14.05 

120 

1.000 

-0.5000 

.577 

.0012 

.0325 

15.22 

130 

1.000 

-0.6428 

.466 

. 0010 

.0280 

16.24 

140 

1.000 

-0.7660. 

... .364 

. 0008 

.0230 

17.08 

150 

1.000 

-0.8660 

.268 

.0006 

.0178 

17.95 

160 

1.000 

-0.9397 

.176 

.0004 

.0122 

18.73 

170 

1.000 

-0.9848 

.088 

..0002 

.0064 

19.65 

180 

1.000 

-1.000 

0 

0 ' 

0 

19.97 

















